Interplay of Coulomb interaction and spin-orbit effects in multi-level quantum dots 
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We study electron transport through a multi-level quantum dot with Rashba spin-orbit interaction 
in the presence of local Coulomb repulsion. We focus on the parameter regime in which the level 
spacing is larger than the level broadening. Motivated by recent experiments, we compute the level 
splitting induced by the spin-orbit interaction at finite Zeeman fields B, which provides a measure 
of the renormalized spin-orbit energy. This level splitting is responsible for the suppression of the 
Kondo ridges at finite B characteristic for the multi-level structure. In addition, the dependence of 
renormalized g-factors on the relative orientation of the applied B field and the spin-orbit direction 
following two different protocols used in experiments is investigated. 
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I. INTRODUCTION 

In linear response transport through quantum dots at 
low temperatures, a two-fold Kramers degeneracy leads 
to the spin Kondo effect in presence of a sufSciently 
strong local interaction U as compared to the level-lead 
hybridization F. In the Kondo regime charge fluctua- 
tions of the dot are suppressed and the physics is dom- 
inated by spin fluctuations. Varying the level positions 
by an external gate voltage Vq, characteristic conduc- 
tance plateaus, so-called Kondo ridges, of width U ap- 
pear around odd (average) electron fillings.l^M^ Breaking 
the two-fold Kramers-degeneracy by a local Zeeman field 
of amplitude B destroys the Kondo ridge and the conduc- 
tance plateau is split into two Lorentzian resonances of 
width oc F along the Vc-axis. In contrast, spin-orbit in- 
teraction (SOI), although breaking spin-rotational sym- 
metry by designating a certain (spin) direction, does not 
destroy the Kondo effect .^^It^ In the presence of SOI spin 
is no longer a good quantum number but a Kramers dou- 
blet remains as time-reversal symmetry is conserved. 

In multi-level dots with initially (at B = 0) well sepa- 
rated levels increasing B might lead to energetically de- 
generate states (level crossings) resulting from different 
orbitals. If one is a spin-up and one a spin-down state 
and the gate voltage is tuned such that an electron fluc- 
tuates between these states one might expect the emer- 
gence of a spin Kondo effect at finite magnetic fields.'^i^tl^ 
If the orbital quantum number is conserved in the leads 
in such systems additional orbital Kondo effect^i^ and 
combinations of spin and orbital Kondo effectJ^^ may 
appear. Here we consider a setup where the dot orbital 
quantum number does not arise in the leads and we thus 
concentrate on the spin Kondo effect. In contrast to the 
standard B — Kondo effect the one appearing at finite 
B is not protected by time-reversal symmetry and can 
be suppressed in presence of a finite SOI.^ 

We here study the dot setup sketched in Fig. [T] as 
a minimal model for a multi-level system. It consists 




FIG. 1. (Color online) The considered setup consists of two 
parallel quantum dots with energies ei/2 = Vg ±5 coupled by 
a hopping amplitude t and a SOI of strength a. The levels 
are split by an external Zeeman field B. The local Coulomb 
interaction is U and the interaction between electrons on the 
two dots is U' . The system is coupled to noninteracting leads 
by hopping amplitudes tpj with P — L,R and j — 1,2. 

of a tight-binding model with two lattice sites 1 and 
2 coupled by the electron hopping of amplitude t and 
connected to two semi-infinite noninteracting leads via 
tunnel couplings of strength t/^^ with /3 = L, R and 
j = 1,2. The on-site energies of the two levels are given 
by ei/2 = ± 5. The Rashba SOI identifies the z- 
direction of the spin space and is modeled as an imagi- 
nary electron hoppi ng w ith spin-dependent sign between 
the two lattice sites.'^EsHU] choose the parameters in 
such a way as to deal with two well separated sets of spin 
degenerate states, that is the molecular regime. We here 
exclusively consider the coupling of a magnetic field to 
the spin degree of freedom (Zeeman term) and neglect 
its effect on the orbital motion. The Zeeman field can 
be decomposed in a parallel and an orthogonal compo- 
nent with respect to the z-direction. The local Coulomb 
interaction (charging energy) is modeled as an on-site U 
as well as a nearest-neighbor interaction U' , and treated 
within an approximate static functional-renormalization 
group (fRG) approach.!^ 

Motivated by recent experiments on InAs quantum 
dots in the Kondo regimCjI^SiMl -^q determine the level 
splitting induced by the SOI at finite B fields. It pro- 
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vides a measure for the renormalized SOI energy. The 
dependence of this spin-orbit energy on the relative ori- 
entation of th e Zeema n field and the spin-orbit direction 
was measured !22l2ll2lHlil In the interpretation of the data 
Kondo correlations were ignored. We here study how 
those modify the angular dependence of the level split- 
ting. In addition, we investigate how the two-particle 
interaction affects the orientation dependence of effec- 
tive g-factors, as extracted experimentally following two 
distinct protocols. In the first the ga te- voltage depen- 
dence of the linear conduct ance^^EMU jg ugg(j extract 
ffCond, in the second bias spectroscopy is employed lead- 
ing to ffLeveis-^^^ '^^ We wiU show that the SOI en- 
ergy and (7Levcis both have an overall amplitude which 
is renormalized by the two-particle interaction. In addi- 
tion, the functional dependence on the angle between the 
SOI and Zeeman field direction is modified. The g-factor 
5Cond remains unaffected. The renormalization due to 
the two-particle interaction is seen to be competing with 
asymmetry effects. 

The paper is organized as follows. In the next section, 
we introduce our minimal multi-level dot model and re- 
view the basic concepts of the approximate fRG treat- 
ment of the Coulomb interaction. In Sec. IIIII we discuss 
our results. We first assess the potential of our model 
in connection with the fRG to describe the experimen- 
tally observed effects in Sec. |III A[ In Sec. |IIIB| we de- 
termine the angular dependence of the SOI-induced level 
splitting at finite Zeeman field for the simplest possible 
model and provide an intuitive physical picture for the in- 
terpretation of the finite-bias spectroscopy. In Sec. |III C| 
we discuss the extraction of effective g-factors from the 
gate-voltage dependence of the linear conductance and 
the bias spectroscopy as well as their respective angu- 
lar dependencies. Sec. |IIID| deals with an extension of 
the simplest models and examines a more realistic asym- 
metric set of paramet ers, ina king close contact to recent 
experimental data.'^^'^llMlIl Finally, we conclude with a 
short summary. 



II. MODEL AND METHOD 
A. Multi-Level Quantum Dot 

The considered minimal multi-level quantum dot 
model is realized by two spin-degenerate levels (at B = 0) 
with the possibility of electron hoppings between these 
levels as sketched in Fig. [T] The Hamiltonian of the iso- 
lated dot contains several terms 
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The free part 
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with dj ^ being the creation operator of an electron of 
spin a =t,i on the dot site j = 1,2 (Wannier states), 
contains the conventional hopping i > 0, and the on-site 
energies ei/2 = Vg±(5 which can be tuned by an external 
gate voltage Vq . The difference of the on-site energies is 
parametrized by the level splitting 26. The effect of a 
Rashba SOI resulting from spatial confinement is taken 
into account by an i maginary hopping amplitude of spin- 
dependent sign.lSlifiJtiil The Rashba hopping term with 
amplitude a > reads 
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with the third Pauli matrix cr^. This choice corresponds 
to a confinement in y-direction if one starts from a one- 
dimensional system in x-direction.'^ We note in passing 
that other SOI terms similar to i?soi with iaz — >■ 
(Rashba SOI from confinement in z-direction) or ia^ 
ay (Dresselhaus SOI)^* can be included in the model but 
will be omitted for simplicity. The SOI breaks the spin- 
rotational invariance and the Zeeman field can be de- 
composed in a component parallel to the SOI (that is 
in z-direction) and one perpendicular to it.l^ We here 
choose the x-direction such that the (local) Zeeman term 
reads 



a.cr' i = l,2 



(2) 



For (f) — ±7r/2 the SOI and the B-field are (anti-)parallel. 
In this case physical quantities are similar to the a = 



cas d ' I if t is replaced by the effective hopping t^s = 
^/W+~a^. In particular, the Kondo ridges at finite B are 



preserved (see Sec. Ill A for a more detailed discussion) 



The local Coulomb interaction is included by 
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for the local U > and nearest-neighbor U' > interac- 
tions respectively, with rij^cr = dj crdj,(T and nj ~ "-j,cr- 
In principle different local interactions Ui , U2 on the 
two sites can be included but we focus on the case 
U = Ui = 1/2- By subtracting 1/2 from rij^a- in the 
definition of TJjnt the point Vq = corresponds to half- 
filling (even in the presence of Coulomb repulsion) of a 
symmetric serial dot, which will be the main geometry 
under consideration in the following. 

Finally, the dot Hamiltonian is supplemented by a term 
describing two semi-infinite noninteracting leads, which 
we model as one-dimensional tight-binding chains 
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c/3 j> + H.c. , (3) 
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with lead operators c^^^ ^ and equal band width 4r. In 
the following we choose r = 1 as the energy unit. We note 
in passing that it is possible to include Rashba SOI terms 
in the leads. For our case of vanishing magnetic field in 
the leads the effects of such terms enter our calculation 
only in terms of an effective hopping.!^ The consequences 
for the Kondo temperature were discussed in Refs. [5], 
[5^ and [3^. The dot-lead couplings are given by the 
tunnel Hamiltonian 



(4) 



Si=L,Rj = \,2 (J 



with tunnel barriers set by 



For simplicity we here 



consider only real tp^^ in the so-called wide-band limit 
(see e.g. Ref. [22]) in which the tunnel barriers only enter 
in combination with the local lead density of states eval- 
uated at the chemical potential. For our setup we find 
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B. Functional RG 



We briefly review the applied approximation scheme 
which is based on the fRCl^S Recent applications to sys- 
tems with SOI inc lude h omogeneous quantum wircs^'' 23 
and quantum dots.l^^l^ 

Starting point of the fRG scheme is the bare {U = 
U' = 0) propagator Qq of the double dot. The leads are 
projected onto the dot sites and enter via the hybridiza- 

tions Tj = S/sT/jj and 7 = Trpieads I]/3 = 



Z)^ ^/9, 1^/3, 2-'^ In the basis 



{|l,t),|l,i),|2,t),|2,i)} 



(5) 



of single-particle dot states the inverse of the propagator 
in Matsubara frequency space reads 
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with Tj(uj) — Tj sgn(a;) and 7(0;) — 7Sgn(ti;). Within the 
fRG Qo is replaced by a modified propagator Qq which 
suppresses low-energy degrees of freedom below a sharp 
Matsubara frequency cutoff A: 

= 6(1^1 -A)goM • 

The cutoff A is sent from 00 down to 0, at which the 
cutoff-free problem is restored. Inserting in the gen- 
erating functional of the one-particle irreducible vertex 
functions, an infinite hierarchy of coupled differential 
equations is obtained by differentiating the generating 
functional with respect to A and expanding it in powers 
of the external fields. Practical implementations require 
a truncation of the flow equation hierarchy. 

Following Ref. [T7], we restrict the present analysis to 
the first order in the hierarchy and only consider the 
flow of the single-particle vertex, that is the self-energy 
S^. Within this truncation T,^ is frequency independent 
leading to a static approximation. It already captures 
the relevant Kondo physics present in the system and 
allows for a qualitative description of equilibrium prop- 
erties such as the linear conductance or the dot occupa- 
tion with minor numerical effort. In the case of a sin- 
gle (spin-degenerate) dot level comparing with numerical 
renormalization group data and Bethe ansatz result^^H 
shows, that the fRG is rehable ff U/T (with F = Fl -FFfl) 
is not too large. In particular, the results for the zero- 
temperature linear conductance and the renormalized ef- 



fective single-particle level show Kondo physics (Kondo 
ridges and pinning of the level energy at the Fermi en- 
ergy). The Zeeman field necessary to suppress the con- 
ductance at half filling to half its value can be used to de- 
fine a Kondo scale Tk- Within first order fRG it is given 
by Tk ^ exp [— ?7/(7rF)] compared to the exact resullP 
Tk ^ exp [— 7rJ7/(8r)]. First order fRG was also shown 
to produce reliable results for multi-level dots as long as 
the same constraint on the ratio of the local Coulomb in- 
teraction and the hybridization as above is fulfilled and 
the number of degenerate single-particle levels does not 
become too large.'^S The quantitative accuracy can be 
improved by including the flow of the st atic par t of the 
two-particle vertex (effective interaction) p2 | 39 |40] ■yyj-^iQ]-^ jg 
beyond the scope of the present analysis. For an in depth 
discussion concerning the range of validity of these ap- 
proximations see Ref. [22]. 

The flow equation for the self-energy read^SSl 

cj=±A b,b' 

where the indices a,a',b,b' label the quantum numbers 
Ta',b':a,b IS thc auti-symmctrized two-particle ver- 
tex, and the interacting Green function Q is determined 
by the Dyson equation 

g^(jt^)= [go-i(H-S^]"' . (8) 



The initial condition for Aq — )■ oo is S"^" = Op2l In the 
lowest-order scheme the two-particle vertex is the bare 
anti-symmetrized interaction and reads: 

S. 



Dynamical contributions to are generated only at 
higher orders. As the latter are important for the con- 
ductance at finite temperatures T > 0, the present ap- 
proximation scheme is restricted to T = 0. The correct 
temperature dependence of the (single-dot) Kondo ridge 
is only captured if the flow of a frequency dependent two- 
particle vertex — leading to a flowing frequency dependent 
self-energy — is kept! " I Within our approximation, the 
matrix elements = T,a',a of the self-energy at the 

end of the flow can be interpreted as interaction-induced 
renormalizations to the noninteracting model parameters 
such as the SOI and the on-site energies .'-^-^ Furthermore 
entirely new matrix elements will be generated if per- 
mitted by symmetry. The full propagator including in- 
teraction effects is determined via the Dyson equation 
from which various observables can be computed. 1^ 
While the renormalized effective single-particle energy 
levels as a function of the bare parameters will be dis- 
cussed in detail in Sec. |IIIB[ 
the linear conductance. At T = 
tions vanish and the Kubo formula for the spin-resolved 
conductance assumes a generalized Landauer-Biittiker 
form^ 



we here concentrate on 
current-vertex correc- 



G(J,<t' 



with the effective transmission T^^o-' (0) evaluated at the 
chemical potential. For the present setup the transmis- 
sion is given by the matrix elements of the full propagator 
leading tcp2l 



III. RESULTS 



A. Linear conductance 



The linear (and finite bias) transport characteristics 
allow to access the physics of quantum dots, and are of 
of particular interest in view o f the use o f dot setups as 
information processing devices .'^^'^^'^Ilini jj^ experiments 
large ranges of applied gate voltages and Zeeman fields as 
well as different orientations of the field can be analyzed. 
The fRG was shown to capture the effects of the two- 
particle in teract ion on the linear conductance for gen- 
eral modelJ^^'^U as well as for models specifically tailored 
to describe experimental setups.'^ We here compute the 
conductance for various parameters of our minimal model 
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FIG. 2. (Color online) Schematic representation of the non- 
interacting level structure of the isolated dot as relevant in 
Sec. HITBl 
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FIG. 3. (Color online) Qualitative reproduction of the con- 
duction region of Fig. Ic in Ref. [23j. The chosen system 
parameters are: t = a = 1, 5 = 0.9, <j) = 0.467r, U = U' = 1.5, 
tL,i = -0.5, tL,2 = 0.45, tR,i = 0.4, tR,2 = -0.3.ES1 



to study multi-level dots with sizeable SOI as realized in 
experiments. As using experimental parameters without 
a well established microscopic model is difficult, we will 
choose our parameters in the following in such a way as to 
facilitate the discussion. As a guide for the magnitude of 
the single particle parameters we use Ref. [S^ . The non- 
interacting level structure of the isolated dot is schemati- 
cally shown in Fig.[2]for ci = 62 = Vq = and can be uti- 
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lized to obtain a rough picture of the linear conductance. 
The approximate fRG conserves particle-hole symmetry 
translating into a symmetric linear conductance with re- 
spect to the gate voltage transformation Vq — ^ — VgI^^I 
For vanishing Zeeman field _B = there are two well sep- 
arated spin degenerate energy levels at e = ±foff- When 
those levels cross the chemical potential of the leads, the 
Kondo effect will lead to the characteristic conductance 
plateaus of height Guax = 2e^//i (for a totally symmetric 
serial dot). The width of the plateaus is determined by 
the local Coulomb interaction. In presence of a finite field 
B parallel to the spin-orbit direction additional spin de- 
generate levels occur at Vq =0, = ±toff (dashed lines 
in Fig. [2]) giving rise to Kondo correlations. Similarly to 
the symmetry in the gate voltage, the conductance is also 
invariant under B ^ —B. This symmetry also holds in 
the following more general cases. If the SU(2) spin sym- 
metry is broken by SOI (defining the z-direction in the 
spin space, see Eqs. ([T]) and ([2|) this finite- i? degeneracy 
can be lifted via a Zeeman field component perpendicu- 
lar to the SOI and the resulting anti-crossing (full lines 
in Fig. 2]) suppresses the Kondo effect on an exponen- 
tial scale. For different coupling strengths to the leads 
the maximum conductance on the Kondo plateau is re- 
duced by a factor which depends on the dot parameters. 
Including additional asymmetries in the on-site energies 
ei — £2 = 2 (5 ^ gives rise to finite- field Kondo ridges 
bent with respect to the Vg axis.'^^'^ Aside from the de- 
velopment of the Kondo effect due to the local Coulomb 
interaction, the nearest-neighbor interaction renormal- 
izes the position of conductance resonances. 

With this simplified multi-level quantum dot model we 
can provide a good qualitative description of the vari- 
ous parameter regimes of recent experiments. To exem- 
plify this in Fig. |3] we show an asymmetrically coupled 
dot with different on-site energies in a Zeeman field with 
a perpendicular component to the spin-orbit direction. 
The results feature a reduced conductance and a slanted 
resonance with a pronounced suppression due to the anti- 
crossing of the levels, both characteristic for the asym- 
metric setup. For the chosen parameters (see the caption) 
the B- and Vq- dependence of G strongly resembles the 
experimental data shown in Fig. Ic of Ref . [53] . 

In the following we use the minimal multi-level model 
to compute the gap in the single-particle spectrum 
(see Fig. [2]) for different relative orientations of the 
SOI and the Zeeman field as recently investigated 
experimentally.'2^l2iII3 A similar level splitting for small 
Zeeman fields close to the B ~ Kondo ridges allows 
to determine the effective factor (7lgvc1s- Alternatively, 
5Cond can be extracted from the Coulomb blockade peak 
splitting of the conductance. For both we will compare 
our theoretical results to the experimental ones. In or- 
der to disentangle the effects under consideration from 
asymmetry effects (possibly enhanced or suppressed by 
the interaction) we will focus on a serial double quantum 
dot geometry. Detailed results for the conductance of 
this model were presented in Ref. 1171 
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FIG. 4. (Color online) Bare spin-orbit parameter a depen- 
dence of the generated effective single-particle parameters Bs 
(upper (blue) curves) and j3 (lower (red) curves) for differ- 
ent orientations of the bare magnetic field (full lines </!> = ^, 
dashed lines tf) — ^, dotted lines (j) — The remaining pa- 
rameters are: t = 1, S = 0, B = 1.2, U = U' = 1, tcoup = 0.4. 



B. Spin orbit energy Esoi 

Recent experiments on InAs devices, where SOI is 
relevant, showed Kondo- like fe atures in the linear con- 
ductance at finite Zeeman field.'^^Ml Cotunneling spec- 
troscopy allows to resolve the dependence of the involved 
energy levels on the B field and its orientation relative 
to the SOI direct ion This analysis shows that any fi- 
nite orthogonal component of the Zeeman fields lifts the 
degeneracy of the two states responsible for the Kondo 
plateau. For the generic situation an anti-crossing of 
these two states is observed, varying periodically with 
the orientation of the Zeeman field. The minimal size of 
the gap as a function of the magnitude of the Zeeman 
field defines the energy E'soi (see Fig. [2| . For the theo- 
retical description we consider the symmetrically coupled 
serial double-dot geometry at vanishing gate voltage and 
level splitting - i.e. ei = £2 = 0, tL.i = i_R,2 = ^Coup and 
tL,2 = = in Fig. [Ij This implies Fi = r2 = F. At 
any (fixed) cutoff value during the RG flow the single- 
particle part of our system (omitting the lead terms of 
Eqs. ([3| and Q) can be described by the following Hamil- 
tonian in the basis of Eq. ([5]): 



B., - iBs 
—t — ia 
il3 



B^ + iBs 
-B, 

—t -I- ia 



—t + ia 
-il3 
B, 
B^ + iBs 



-il3 
—t — ia 
B^ - iBs 
-B, 



with all matrix elements depending on the cutoff A. The 
initial conditions are given by the bare values of the non- 
interacting system, with Bz = B sin (f> and B^ = B cos (f). 
The Hamiltonian (matrix) h contains the parameters Bg 
and /3 which are zero initially but are generated by the 
two-particle interaction during the RG flow. The new 
parameter Bs is a Zeeman field perpendicular to both 
the applied Zeeman field and the SOI direction, with op- 
posite orientation on the two dot sites, while the spin 
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FIG. 5. (Color online) Bare Zeeman field dependence of the 
effective single-particle energy spectrum (full lines) as well as 
of the corresponding noninteracting ones (dashed lines). As 
a guide to the eye the chemical potential and the crossing 
position of the renormalized levels are given as straight black 
lines. The parameters are: Vg = 1.25, t = 1, a = 0.6, (j) = 



0.257r, U = U' = 1, tcou 



0.4. 



flip hopping /3 is a Dresselhaus SOI term.'i^ The ap- 
pearance of an effective Zeeman field induced by a finite 
Coulomb interaction in presence of a broken spin sym- 
metr y is d iscussed for quantum dots with ferromagnetic 
leadP^I^ a nd has been observed recently in systems in- 
volving SOI.ESm We find a non-monotonic dependence 
of the effective IBs = B^^°, (3 = for varying the 

initial value of the SOI - given by a - and initial relative 
orientation of the Zeeman field and the SOI parametrized 
by the angle cj) as shown in Fig. [4j 

At the end of the flow the effective parameters in- 
clude the renormalization of the initial values due to the 
Coulomb interaction. Here we focus on the angular de- 
pendence. For the considered setup the flow equations 
assume the convenient analytical form in terms of vec- 
tors B, = {B^,B,,BsV and 4 = (a, /3, tf 



B, 



U 



U' 



ttD{A) 
where we introduced 



2 ytc ■ Bcj tc + f+Bc 
2(t,-B,^B, + f_t, 



(9) 



/ = |A| + r, /± = /2±(| 



Br 



D{A) det(i/ ~h)^fl + 2f\Q^ + 4 fe • Be 



As a consequence, Bs is generated only for both finite 
initial B^ and a (compare to Fig. [4]) . For the hopping 
/3 to be generated, a non-vanishing Bz is additionally 
required. 

Diagonalizing the Hamiltonian h yields the eigenval- 
ues (the symmetric leads contribute a constant imaginary 



FIG. 6. (Color online) Renormalization-induced (j) depen- 
dence of the effective parameters. Here we compensated for 
the bare angular dependence of the Zeeman fields by divid- 
ing by the corresponding trigonometric function. The initial 
parameters are: t — 1, a — 0.6, B — 1.2, U — U' = 1, 
tcoup = 0.4. 
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FIG. 7. (Color online) Spin-orbit energy versus orientation 
of the Zeeman field for several values of the interaction U = 
)7'(full lines). The dashed lines show the spin-orbit energies 
normalized to the angular dependence of the noninteracting 
case. The initial conditions are given by f = 1, a — 0.6, tL = 
tR = 0.3. 



part which is omitted here, compare Sec. HID) 



|i3J2±2W|BJ2|t 



They can be interpreted as the single-particle levels of 
a corresponding noninteracting system. This interpreta- 
tion has already been applied successfully to the problem 
of phase lapses in multi-level quantum dots.'^" Figurejs] 
shows data for finite gate voltages for which pronounced 
renormalization effects occur if one of the involved levels 
crosses the leads' chemical potential. For this situation 
the other energy levels are shifted upwards due to charg- 
ing effects. This is not observed at Vq = as no level 
crosses the chemical potential, but the overall shift of the 
effective levels with respect to the noninteracting ones ap- 
parent in Fig. [5] remains and will be of importance in the 
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following. From the two intermediate levels we determine 
the spin-orbit energy 

Esoi = 



mins j2y|reP + |i3eP-2y|4P|roP-(re-4) j . 

The bare initial value -Esoi = 2 a cos (f> — 2 a± is ob- 
tained for -B^i„ = -|-a^|, with q!|| = a sin^. Due to the 
complicated non-linear structure of the flow equations an 
analytic expression of the renormalized spin-orbit energy 
in terms of the bare Zeeman field can not be obtained. 
The numerical solution of Eq. ([t]) or ([9| shows that the 
renormalized effective parameters acquire a non-trivial 
angular dependence, as seen in Fig.|6]The results depend 
only quantitatively on the details of the Coulomb inter- 
action and the inter-dot hopping t, as long as C/^'^ /T and 
t/r are sufficiently large. We here focus on U — U' and 
t = 1. The renormalized spin-orbit energies are shown 
as solid lines in Fig. [7] and follow the general form of the 
bare case calculated above: We find a maximum of the 
level splitting for the perpendicular orientation (f> — 
and a monotonous decrease to i?soi = when increasing 
(j) towards the parallel configuration at (j) — tt /2. In a de- 
tailed examination for more values of the interaction the 
maximum is seen to increase linearly with U . The simi- 
larity of the interacting and noninteracting curves might 
lead to the expectation that the maximum is given by 
twice the renormalized SOI parameter 2 a, but this is 
not the case. The dashed lines in Fig. [7] show the renor- 
malized spin-orbit energies divided by the angular depen- 
dence of the bare case, i.e. cos 0. The strong deviations 
from this bare dependence close to the parallel orien- 
tation leads to the important result that the functional 
dependence of the SOI energy on (j) is strongly affected by 
the two-particle interaction. This can be traced back to 
the more pronounced renormalization effects of the bare 
parameters around the parallel configuration as is seen 
in Fig. |6j while the parameters are mostly unaffected 
for (j) ^ O.Stt. This behavior can be intuitively under- 
stood in the following way: for the parallel configuration 
a sufficiently strong Coulomb interaction gives rise to the 
finite-B Kondo effect related to the level crossing. With 
the vanishing Kondo effect in presence of a finite orthog- 
onal S-field component relative to the SOI the Coulomb 
interaction effects appear suppressed as well. We finally 
note that our results agree qualitatively with the experi- 
mental results of Refs. [23] and [5S]. 

C. Effective g-factors 

Experimentally the (/)-dependence of effective g-factors 
is studied as well. In this section we will model two dif- 
ferent protocols used for their extraction. In the first one, 
following the experiments of Refs. [53], [M], and [35] the 
Coulomb blockade peak splitting around a, B = Kondo 
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FIG. 8. (Color online) Determination of the g-factor gcond 
from the gate-voltage dependence of the Coulomb blockade 
peaks in the linear conductance. Bare parameters: t — 1, 
a = 1, (j> ~ O.Stt, U — U' = 0.5, tcoup = 0.3. The inset shows 
a linear fit to the extracted gate voltage difference A. 

resonance is extracted from the linear conductance defin- 
ing gcond- As the fRG reliably reproduces the linear con- 
ductance this procedure is easily adopted. In the second 
protocol bias spectroscopy is used to measure the level 
splitting in the vicinity of a Kondo resonance, determin- 
ing gLcvols- 



Effective g-factors: gcond 

As seen in linear conductance measurements for B > 
resonance peaks of maximal height /h and width ^ F 
develop out of Kondo plateaus, corresponding to the fill- 
ing of a dot state. From a linear fit to the splitting of 
these Coulomb blockade peaks at small to intermediate 
Zeeman fields we determine gcond by identifying it with 
the slope of the fit as shown in Fig. |8| see i.e. the experi- 
ments reported in Ref. [S^, Fig. Ic. Due to the presence 
of the B — Kondo ridge, there is an offset before the 
linear behavior sets in. The full angular dependence of 
5Cond in the interacting system is shown in Fig. |9| and 
for comparison also the noninteracting result is displayed. 
No significant effects of the two-particle interaction are 
observed. For all values of the interaction, gcond exhibits 
an S-shaped dependence on the relative orientation of 
the Zeeman field and the SOI with a maximum for the 
parallel configuration and a finite minimum at (j) — 0. 



Effective g-factors: gLcvds 

In order to use the fRG to model the effective g-factor 
^Levels we will make use of the effective-level interpreta- 
tion, according to the discussion in the previous section.l^Sl 
Introducing a finite gate voltage as additional parame- 
ter, an explicit analytic form of the flow equation as for 
the spin-orbit energy (see Eq. ([9])) is much more difficult 
to obtain. Thus we will solve the general fiow equation 



FIG. 9. (Color online) Zeeman field orientation dependence 
of (?Cond as extracted from the gate-voltage dependence (see 
Fig. |8| for different values of the Coulomb interaction U = 
U' = 0, 0.25, 0.5, 1 and parameters as in Fig. [s] 



Eq. ([t]) including a numerical inversion of the matrix on 
the right hand side of Eq. ([8| and extract the low field 
splitting from the eigenvalues of the resulting effective 
Hamiltonian. As for gcond, we follow the experimental 
procedur^^, perform a linear fit to the computed split- 
ting and identify the slope as ^Levels- In principle we can 
extract a gate-voltage dependence for this quantity but 
here we choose the gate voltage such that the linear con- 
ductance around the B — Kondo resonances is maxi- 
mal. We note that in the single-impurity Anderson model 
^Levels Can bc related to the magnetic susceptibility.l^ 

The noninteracting Hamiltonian with ei = £2 = Vg 
yields the eigenvalues 



Vg± 



'II 



from which the bare gLcvois is identified 



|Ae±| 



'II 



cfi' 



cvcls 



B 



for B < f 



The angular dependence for the symmetric case at fixed 
Vq is shown in Fig.[TO]for different values of the Coulomb 
interaction. The general form is again S-shaped as for 
ffCond- In contrast to gcond, for ^lgvoIs we find similar in- 
teraction effects as for -Esoi- Again the renormalization 
of the maximum value of 5lcvo1s is a prominent effect. 
Comparing the 0-dependence of the interacting curves 
with the noninteracting ones we find deviations which 
are much less pronounced compared to the ones of the 
spin-orbit energy (compare dashed lines of Figs. 10 and 
[7|. Even though the qualitative 0-dependence is similar 
to the one of gcond, the strong renormalization of the am- 
plitude implies that (^Levels 7^ 5Cond- This makes it neces- 
sary to clearly distinguish between these two quantities 
as well as other contextually similar definitions using dif- 



FIG. 10. (Color online) Angular dependence of ^Levels as 
extracted from the level splitting (solid lines) for the same 
parameters as in Fig. |8] and different values of the Coulomb 
interaction U = U' — 0,0.25,0.5,1. In addition, the angular 
dependence normalized to the noninteracting case is shown 
(dashed lines). 




FIG. 11. (Color online) Angular dependence of the spin-orbit 
energy for the parameters of Fig. |3] (full lines). For reference, 
the angular dependence of the noninteracting case is shown 
(dashed lines). The inset shows a zoom in of the (j) > OAir 
range, note the logarithmic scale on the y-axis. 



ferent computation methods or extraction prof o cols .l^^^^^ 



D. Asymmetry effects 

Experimental setups are more accurately modeled if we 
consider the more general situation of four (asymmetric) 
couplings and different on-site potentials as depicted in 
Fig. 1. Conceptually and computationally this is eas- 
ily implemented within the fRG approach but due to the 
multitude of parameters an in depth analysis is beyond 
the scope of this work. Instead we will use the experimen- 
tally motivated parameters of Fig. 3 and calculate the 
spin-orbit energy and effective g-factors for this asym- 
metric setup to illustrate asymmetry effects. 
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Spin orbit energy Ssoi 

While the calculation of the g-factors remains unaf- 
fected by the setup, the protocol for the spin-orbit en- 
ergy has to be slightly expanded. In the experiments the 
parameters were tuned in such a way that the dot was at 
half filling to insure maximum degeneracy. In Sec. |IIIB| 
this half filling condition was guaranteed for Vq = due 
to particle-hole symmetry. In the considered asymmet- 
ric setup and for finite two-particle interaction this does 
not hold necessarily and we will fix the gate voltage and 
Zeeman field amplitude so that on average two electrons 
occupy the dot and the considered level splitting Esoi 
is minimal. For the noninteracting case the dot is still 
half filled at Vq = and the spin-orbit energy can be 
computed after diagonalizing Eq. Q for a; = 0+. As 
Fi 7^ F2 and 7 7^ the lead self-energy contribution 
affects the real (and the imaginary) part of the eigenval- 
ues of the dot system. In Sec. |IIIB| the lead self-energy 
contribution was proportional to unity and thus the real 
parts of the eigenenergies of the isolated dot remained 
unaffected. This renormalization by the leads has a dras- 
tic effect on the spin-orbit energy as the dashed lines in 
Fig. fill show. The overall structure remains similar to 



2.5 



the serial symmetric case of Sec. |IIIB with a reduced 
maximum (compared to 2 a in Sec. Ill B[ ) at (j) — and 
a decrease towards larger (j). The spin-orbit energy then 
tends towards for « 0.457r and remains very small 
up until (f) = O.Stt. The picture for the interacting dot is 
similar (full line in Fig. 11) but the two-particle interac- 



tion decreases the range for which the spin-orbit energy 
is small compared to the U = case. In the theory plots 
this lead renormalization effect occurs most visibly for 
large (fi where the spin-orbit energy is small. We refrain 
from comparing this effect to the experiment as it might 
be masked by finite temperature effects or the resolution 
in the bias spectroscopy. It is obvious that this effect de- 
pends strongly on the detailed geometry of the quantum 
dot. A more thorough analysis is needed to determine 
the importance of this result for experiments. 



Effective g-f actors 

For gcond we find several interesting features compared 
to the symmetric case in the previous section. The most 
obvious effect seen in Fig.[T2]is a different strength of the 
effective y-factor depending on which B = level degen- 
eracy (U = 0) or Kondo resonance {U ^ 0) is selected for 
measurement (Vq = g*). This is already seen in the 
noninteracting curves depicted as dashed lines in Fig. [12] 
and is only weakly affected by the interaction (full lines). 
An interesting effect is observed if one considers the an- 
gular dependence of the noninteractin g case . Only g^^^-^^ 
follows the general form found in Sec 



IIIC 



while 

remains nearly constant over the whole </> range (slight 
deviations might be attributed to numerics). If the in- 
teraction is turned on this atypical behavior is not found 
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FIG. 12. (Color online) Angular dependence of g^^^^ for the 
parameters of Fig. [s] (full lines). For reference the dependence 
of the noninteracting case is shown by the dashed lines. 




FIG. 13. (Color online) Angular dependence of gi^^^^i^ for the 
parameters of Fig. [s] (full lines). Normalization with respect 
to the bare angular dependence of the noninteracting case is 
shown by the dashed lines. 



and both curves follow the general form with only slightly 
renormalized amplitude compared to the free case. 

Computin g ffLcvc is we find that the angular dependence 
found in Sec. Ill C is qualitatively preserved (full lines in 
Fig. 13). The two-particle interaction renormalization 



effect on the detailed functional form of the angular de- 
pendence is suppressed in comparison to the symmetric 
setup (dashed lines in Fig. 13). As for gcond a strong de- 
pendence on the asymmetry is seen if we consider both 
B = Kondo ridges, but only in the interacting case. 
As the two sets of states involved in forming the Kondo 
effect are coupled differently to the leads, the renormal- 
ization of parameters close to the resonance is different 
as well. This results in different ^Levels ^^'^ Kondo 
plateaus at positive and negative gate voltages. This 
asymmetry effect present in gcond as well as gLevcis has 
already been observed in experiment^^HMlll where the 
two Kondo ridges can be attributed to different orbitals 
of the device. 
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IV. CONCLUSION 



We studied the influence of the Coulomb interaction 
on the level splitting induced by the spin-orbit interac- 
tion (the spin-orbit energy), and the efl[ective g- factors 
in multi-level quantum dots with SOI. Furthermore, we 
interpret fRG results in terms of effective single-particle 
energy levels to obtain an intuitive physical picture for 
the understanding of finite-bias spectroscopyP^ For a 
basic symmetrically coupled serial model we find that 
two experimentally investigated quantities are affected 
by the local Coulomb interaction. In the case of the 
spin-orbit energy Egoi the effect on the overall ampli- 
tude as well as the dependence on the relative orienta- 
tion between the SOI and the applied Zeeman field is very 
pronounced even for intermediate interaction strengths. 
For the gLcvois-factor calculated from the gate- voltage de- 
pendent effective level splitting - which mimics the ex- 
perimental g-factor extraction from bias spectroscopy - 
we find sizable renormalization of its magnitude while 
the angular dependence is only mildly affected. On the 
other hand, the gcond-fa-ctor extracted from the Coulomb 



blockade peak splitting of the linear conductance appears 
to be almost interaction independent. The considered 
asymmetric parameter set shows qualitatively similar re- 
sults. They nevertheless reveal a complex interplay of 
lead coupling, two-particle interaction and SOI. The pre- 
sented results are of importance for the understanding 
of transport measurements of multi-level quantum dots 
with SOI in presence of an external Zeeman field as re- 
ported in Refs. [23], [H], [IB], and [SI]. While the quali- 
tative behavior of the considered quantities is consistent 
with the experimental data, the observed deviations and 
interaction-dependent amplitudes in the theoretical cal- 
culations provide new directions to investigate in future 
experiments. 
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